Silent Phase Qubit Based on rf-Wave Josephson Junctions 
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We report on design and fabrication of a new type of flux qubit that capitalizes on intrinsic 
properties of submicron YBCO grain boundary junctions. The operating point is protected from the 
fluctuations of the external fields, already on the classical level; the effects of external perturbations 
are absent, to the second or third order, depending on the character of the coupling. We propose 
an experiment to observe quantum tunneling and Rabi oscillations in the qubit. Estimate of the 
decoherence due to fluctuations of the external flux is presented. 
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Over the last few years a series of experiments^* 2 ^ pro- 
vided a conclusive evidence of quantum superposition in 
meso- and macroscopic superconductors. The task at 
hand is scaling up of the system, with two goals in mind: 
to probe how far quantum superposition can be pushed 
into macroscopic world, and to develop an element base 
for quantum computing, which only becomes viable on 
the scale 10-100 qubits. 

It was suggested, that use of high-T c cuprates in qubits 
would dispose of the need to apply fine-tuned exter- 
nal fields to keep it in the operating point, due to the 
time-reversal symmetry breaking in systems with DD 
junctions^. On the level of a few qubits this is not a 
major advantage, compared to the relative difficulty of 
fabrication and threat of extra decoherence from nodal 
quasiparticles and zero-energy states (ZES) in cuprates. 
Therefore the research was concentrated on conventional 
superconductors, where all of the aforementioned suc- 
cesses were achieved. 

Nevertheless, the logical progression of research from 
single qubit to qubit-qubit coupling and further to qubit 
registers is bringing us to the point where intrinsic bista- 
bility of high-T c qubits will become a major advantage. 
In the last few years, we saw the development of reliable 
fabrication of bistable submicron high-T c structures^iS. 
The danger of nodal quasiparticles and ZES seems to 
be overestimatedSiiS. Finally, putting aside any techno- 
logical motivation, it is tempting to probe high-T c com- 
pounds for the same kind of macroscopic quantum coher- 
ence as was observed in low-T c superconductors, just to 
see how and whether these phenomena depend on, e.g., 
symmetry of the ground state of the system. 

The "silent qubit" (for the etimology, see below) can 
be viewed as a SQUID formed out of high-T c film, with 
two DD-grain boundary junctions forming a mesoscopic 
island (Fig. QJi,c). The crystal lattice and d-wave or- 
der parameter orientations on the two sides of the grain 
boundary are chosen in such a way that each junction 
has doubly degenerate ground stat o?' 11 . 

To simplify the analysis, we keep only two harmon- 
ics in the current-phase relations of the DD junctions: 
h = Id sin 4>i ~~ I'd sm 20^, where I C i is the Josephson cur- 
rent, and d>i is the phase difference across the i th junc- 
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FIG. 1: (a) "Silent" qubit. The dashed line shows the loca- 
tion of the grain boundary junction, (b) Suggested scheme for 
observation of Rabi oscillation and measurement of decoher- 
ence. (c) AFM picture of a fabricated qubit with a dc-SQUID 
fabricated around it for readout of the states. 



tion [i = 1,2). This approximation successfully repro- 
duces the observed behaviour of similar devices in classi- 
cal regime^. The flux quantization condition binds the 
phases (pi and 02 to the total flux through the loop $ 
via <j> = (f>i + (j)2 — 27r$/$o, where $ = h/2e is the 
flux quantum. Introducing the superconducting phase of 
the island = (0i — </>2)/2, the free energy of the qubit, 
threaded by an external magnetic flux <fr x , in units of the 
Josephson energy Ej (= E\, where E{ = Id^o/^Tr), is 
given by 
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Here a t = 2I' ci /I ci , 
[(ax +rja 2 ) cos (£] /[(! 
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2-k<& x /<&q. The dimensionless self-inductance of the 
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loop Pl = 27rL/ c i/$o is considered negligible (0l —> 0); 
then <p — > <f> x and the first term in can be dropped. 
The last term in ffl: 
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is zero when the two junctions are identical (i.e. r\ = 1 
and ai = 02). In this case the potential energy min- 
ima are at 9 = ± arccosfd/a^) when > 1 and zero 
otherwise. The ± signs correspond to the states on the 
right and left sides of the potential well (i.e. |±)), respec- 
tively. For a<f, > 1 the potential has two minima, which 
are degenerate at any external flux cj) x . The current in- 
duced in the loop by the external flux does not depend 
on the state of the qubit, which justifies the moniker 
"silent" . The potential profile (JTJ is similar to the one of 
persistent current qubit^, but here we have only one in- 
dependent phase (as /3l — > 0), and the problem becomes 
one-dimensional. The barrier between the potential min- 
ima is flux-dependent: W = cos((j> x /2)(a l p + ai 1 — 2). 

In the general case (ai 7^ 012, T) 7^ 1) the two minima 
are only degenerate when <f> x — 0. Now a state- dependent 
persistent current flows in the lo op even in zero external 
field; in units of I c i, Iq = ±[i]y / al — 1/(1 + vi)a^](a2 — 
ai), where fio = a0=o- 

An intermediate regime (a.\ = a.%, rj ^ 1) is most 
interesting. It takes place when junctions have similar 
current-phase dependences, but different critical currents 
(i.e. widths), and should be expected if the junctions are 
close enough to each other (Fig. QJ. At <f> x — 0, the 
equilibrium value of 9 is the same for both junctions, 
and there is no spontaneous current: Zq =0. At finite 
<j) x , the induced currents differ for the two states of the 
qubit, but the difference is of higher order in <fi x , keeping 
the qubit silent. 

Expanding the free energy to the third order in (j) x , we 
find for the minima: 
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where Ai are explicit, but cumbersome, functions of &o 
and rj. As expected, there is no first order dependence 
on (j> x (A\ = 0). The second order term in @ does 
not depend on the state of the qubit. The first state- 
dependent term in the minimum energy of the system 
is 0(cf> x ) (Fig. Therefore small fluctuations of cf> x do 
not affect the degeneracy of the states. The difference 
between the energy minima grows as the external flux is 
increased until the point at which the potential barrier 
vanishes altogether, and the minimum with higher energy 
disappears (the jumps in Fig. EI> - 

The current in the loop is found from = d^^U^ = 
2A2(j) x ±3A3(f) x . This current generates a state-dependent 
magnetic flux 5(f> = 4> — <fi x = PlI^ ■ Note that the state- 
dependent contribution to the induced flux is 0((f> x ). For 
finite self-inductances /3l, S(f> has to be calculated self- 
consistently. Figure 03 shows the result of such a cal- 
culation. With the parameters chosen, an external flux 
close to 0.2<I>o generates an additional flux (~ 0.005$o) 
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FIG. 2: The minimum of the free energy (a) and the corre- 
sponding values of 6 at the minimum points (b) as a function 
of 4> x , for n = 2 and ot\ — 012 = 3. The solid and dashed lines 
correspond to the | — ) and |+) states respectively. 



through the loop, which is of the same order as the esti- 
mate based on the above expansion. 

Tunneling between the potential minima occurs due 
to the uncertainty relation between the charge Q of the 
island and its superconducting phase 9. The tunneling 
matrix element is approximately given by (h = &b = 1) 

A pa 



t> x )e~^ cw ^ Ej / E ", where ( is a constant of the 



order of 1. 
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The coefficient u) p (4> x ) = 
mined by the frequencies uj^ of small oscillations in the 
right and left potential minima, respectively. This de- 
pendence follows from the expression of A as the matrix 
element between the lowest energy states in the two wells. 
In the case of &o ~ 1 , it is only valid qualitativelj*^, but 
enough for the present analysis. Due to the symmetry 
of the potential profile when evi = 02, the linear de- 
pendence on cj) x cancel and we are left with u> p (<f> x ) = 

^EjE c (ao - a^)(l -«</£), where E c = e 2 /2C is the 
charging energy, C is the effective capacitance of the junc- 
tions, and k is a dimensionless coefficient of O(l). Fluc- 
tuations of 4> x influence W and therefore A. Expand- 
ing the Josephson potential near the origin, we obtain 
the tunnelling barrier W = W ma x - W^ in = (Bo - A ) + 
(B 2 - At) <P 2 X + 0((g), where B = - (r, + 1) (1 - a /4) , 
B 2 = — [(<So — l)/4] [17/(7] + 1)]. Again there is no depen- 
dence on <p x in the lowest order. 

Truncating the Hilbert space of the qubit to the two 
lowest energy states, one can write the effective Hamilto- 
nian of the system as H — ^A(<j) x )a x + he((j) x )cr z , where 
e((f> x ) s» EjA 3 cj) x . All single qubit operations can be real- 
ized by applying controlled flux cf> x . Note that the qubit 
only leaves the operating point when a finite external flux 
is applied. Unlike the earlier qubit designs^ this point 
is protected from external flux fluctuations, already on 
the classical level [cf. Eq. 

Readout of the quantum state can in principle be 
achieved using a dc-SQUID to measure the magnetic flux 
generated by the qubit when it is biased by an external 
magnetic field. In this work the readout SQUID was fab- 
ricated onchip using a high-T c grain boundary SQUID 
(cf. Fig.^:); alternatively it can be made of aluminum. 
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FIG. 3: Self-consistent calculation of the self-generated flux 
for the system of Fig. H with f3 L = 0.01. 

Although coherent oscillations in an all- aluminum per- 
sistent current qubit have previously been observed us- 
ing a similar readout design^, the measured decoherence 
time appeared to be rather short. Subsequent results of 
Il'ichev et al>i^ have demonstrated a longer decoherence 
time in a measurement scheme which couples the qubit 
to a high quality tank circuit. Here, we calculate the de- 
coherence time of our silent qubit in the same setupi&i^ 
(see Fig. [IJ}) . The circuit is assumed to have a resonance 
frequency lot and the damping rate yr- The mutual 
inductance is M — k^/TrpL, k < 1 being the coupling 
coefficient, and L and Lt are the inductance of the qubit 
and tank respectively. 

Without limiting the generality of our approach, we 
can ascribe all the dephasing and dissipation in the 
qubit to its interaction with the fluctuating flux (f> x {t), 
created by the tank current Ir{t)- These fluctuations 
are characterized by correlator K(t,t') — (7r(t), 7r(t')), 
spectral density K(co) = Lo('fTLO T /LT)[coth(Lo/2T) + 
1]/[(lu 2 — lo t ) 2 + y T L0 2 ], and dispersion (if,) = K(t,t) 
= (w T /2iT)coth(w T /2T). At <j) x = 0, the qubit Hamil- 
tonian becomes H — (A/2)a x — X 3 I t <t z — \il T a x , 



where A3 and A2 are coupling coefficients depend- 
ing on the qubit parameters. In the Bloch-Redfield 
approximation^, we calculate the energy relaxation rate 
r = 30A|(/|) 2 (7 T /A 3 )(cj|,/ J L T )coth(A/2T), together 
with a dephasing time of the qubit 7 _1 , where 7 = 
r/2 + 7o , 7o = {16n/3)\ 2 (y T T 3 /L T LO T ), if T < co T , and 
70 = \\Q 't (wj / L 4)[coth 2 {lot /2T) — 1] at temperatures 

T > LOT- 

Using the experimental data of Ref. ^3 {Ic — 0.5 
fj,A, I' c = 0.6 fiA, C « 10 fF), we find E c « 2 GHz, 
io p /2tt w 40 GHz, and A/2tt w 1.6 GHz. For (3 L ~ 0.01, 
which is the value used in Fig. [31 the inductance of the 
loop is of the order of lOpH. To estimate the contribu- 
tions of the cubic and quadratic terms to qubit dephas- 
ing and dissipation, we chose the following parameters: 
T) = 2, a t = a 2 = 2.4, and Ej = 1.66 x 10" 22 J. If 
the tank frequency lot /2k =10 MHz, its quality factor 
Qt = 2000, and the coupling coefficient k ~ 1/33, then 
contribution of quadratic flux fluctuations to dephasing 
rate is small, so that the dephasing time due to qubit cou- 
pling to the tank is y^ 1 ~ 20 ms at temperatures of order 
10 mK, while the contribution of the cubic fluctuations 
to the dephasing and relaxation rates is totally negligi- 
ble. It means that at the operating point the silent qubit 
is practically decoupled from the fluctuations caused by 
the controlling circuits. The dominant source of deco- 
herence is from the nodal quasiparticles at the junction, 
considered in Ref. 0, which may reduce the decoherence 
time to about 1-100 ns. 

In conclusion, a new type of flux qubit, based on spe- 
cific properties of submicron YBCO grain boundary junc- 
tions, is proposed and fabricated. The symmetry of the 
device provides an operating point, which is stable and 
protected against the external field fluctuations. 

We are grateful to A. Golubov, E. Il'ichev, A. Maassen 
van den Brink, Y. Nakamura, and A. Shnirman for stim- 
ulating discussions. 



1 Y. Nakamura, Yu. A. Pashkin, J.S. Tsai, Nature 398, 786 
(1999); J. Friedman et al., Nature 406, 43 (2000); J.M. 
Martinis, S. Nam, and J. Aumentado, Phys. Rev. Lett. 
89, 117901 (2002); Y. Yu, S. Han, X. Chu, and Z. Wang, 
Science 296, 889 (2002); A. J. Berkley et al, Science 300, 
1548 (2003). 

2 I. Chiorescu et al, Science 299, 1869 (2003). 

3 D. Vion et al., Science 296, 889 (2002). 

4 L.B. Ioffe et al, Nature 398, 679 (1999). 

5 A. Zagoskin, |cond-ma t/9903170 A. Blais and A.M. 
Zagoskin, Phys. Rev. A 61, 042308 (2000). 

6 E. Il'ichev et al, Phys. Rev. Lett. 86, 5369 (2001). 

7 A.Ya. Tzalenchuk et al, ASC 2002 (Houston). 

8 T. Lindstrom et al, Phys. Rev. Lett. 90, 117002 (2003). 



9 M.H.S. Amin and A. Yu. Smirnov, preprint 
cond-mat/0304255 

Ya. V. Fominov, A. A. Golubov, and M. Yu. Kupriyanov, 
JETP Lett. 77, 587 (2003). 

1 M.H.S. Amin, A.N. Omelyanchouk, and A.M. Zagoskin, 
Phys. Rev. B 63, 212502 (2001); M.H.S. Amin et al, Phys- 
ica B 318, 162 (2002). 

2 A.Ya. Tzalenchuk et al, Phys. Rev. B 68, 100501(R) 
(2003). 

3 M. Grajcar et al, preprint cond-mat/0303657 

4 E. Il'ichev et al, Phys. Rev. Lett. 91, 097 906 (2003). 

5 AYu. Smirnov, preprint cond-mat/0306004 



